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Abstract 

The mass and isotope dependence of limiting temperatures for hot nuclei are 
investigated. The predicted mass dependence of limiting temperatures is in 
good agreement with data derived from the caloric curve data. The predicted 
isotope distribution of limiting temperatures appears to be a parabolic shape 
and its centroid is not located at the isotope on the /3-stability line(T=0) but 
at neutron-rich side. Our study shows that the mass and isotope dependence 
of limiting temperatures depend on the interaction and the form of surface 
tension and its isopin dependence sensitively. 
PACS numbers: 24.10.-i,25.70.-z,25.70.Gh 



1. INTRODUCTION 

It is well known that for nuclear matter, liquid and gas phases can coexist at tempera- 
ture lower than T c and only gas phase can exist at or above T c . While for a hot nucleus 
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surrounding by vapor, due to the Coulomb instabity [1-4], it may not be stable above a 
temperature which is much lower than the critical temperature T c of infinite nuclear matter, 
this temperature is called the limiting temperature Tu m . Limiting temperatures for nuclei 
along the /^-stability line were studied by using Skyrme interaction [1,2,5-8] and Furnstahl, 
Serot and Tang (FST) model [9]. Their studies showed that the limiting temperature was 
sensitive to the effective nucleon-nucleon interactions. In ref. [10,11] the mass dependence of 
limiting temperatures for hot nuclei was derived by analyzing the existing caloric data. The 
establishment of the beam of rare isotopes will provide us the opportunities to study the 
heavy ion collisions with large isospin asymmetry. The measurements of caloric curves(or 
other possible observables) can possibly explore the limiting temperatures for hot nuclear 
systems with large isospin asymmetry. The limiting temperatures for nuclei away from the 
/3-stability line should depend on the isospin dependent part of the Equation of State(EOS). 
The newly developed Skyrme forces such as the SLy series are designed to study the prop- 
erties of nuclei away from the /3-stability line in addition to nuclei along the /^-stability line. 
It has been shown that these modern Skyrme forces can describe the properties of the nuclei 
away from /3-stability line much better than the old Skyrme forces [12]. Therefore it seems 
to us to be worthwhile to apply the modern parametrization of Skyrme force such as the 
SLy series for studying the limiting temperatures for nuclei both along and away from the 
j3- stability line. 

Additionally, the surface tension is another important factor affecting the limiting tem- 
perature. The surface tension represents a work per unit area needed to create the surface, 
the work is required because the nucleons at the surface are less bounded. Therefore the 
surface tension is related to the EOS of nuclear matter. The symmetry energy term in the 
EOS should also contribute to the surface tension. However, the relationship is not avail- 
able and it is treated independently of the equation of state except the fact that the surface 
tension vanishes as T > T c . In this case, the form of the temperature dependence of the 
surface tension is not unique in literatures and little attention has been paid to the influ- 
ence of the different forms of the surface tension on the limiting temperature. The isospin 
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dependence of the surface tension is even unclear. Therefore it is requisite to investigate the 
effect of the different forms of the surface tension, including with and without the isospin 
dependence(surface symmetry term), on limiting temperature. By comparing the calcula- 
tion results with the experimental data we may obtain the information of the surface tension 
and better parametrization of Skyrme interaction. 

In this work we follow the approach of Ref. [13,14] to study the mass and isotope depen- 
dence of limiting temperatures by using the Skyrme effective force of SLy7 in addition to 
SIII and SkM*. The effect of different forms of the surface tension on the limiting tempera- 
ture will also be investigated. In Sec. II, The equation of state by using Sly7 is studied. In 
Sec. Ill we give the coexistence equations and the calculation results of the mass and isotope 
dependence of limiting temperatures of hot nuclei. Finally, a summary and discussions is 
given in Sec. IV. 

2. THE EQUATION OF STATE 

The single-particle energy reads 

£q = 2m* + Uq + £ ° ou 9 ' p ' ^ ' 

where m* and u q is the effective mass and the single-particle potential energy of species q, 
respectively. The third term is the Coulomb energy. The chemical potential p q of species q 
for nuclear system at temperature T can be written as 

A 3 00 n + 1 A 3 
p q (T, p n , p p ) = u q (p n , p p ) + e C oui5q,p + T In — p q + K{—p q ) n - (2) 

9s n =l 71 9s 

Here At is the effective thermal wavelength of the nucleon, which reads 

,2ti% 2 . i . , 

A T = 3 
m*T 

and bnS are the coefficients of the virial series for ideal Fermi gas. For the most general 
form of a Skyrme-type interaction, u q can be written as 
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Uq = A lP + A 2 p q + A 3 r + A 4 r q + A 5 (a + 2)p a+1 

+A 6 ap a - l (p 2 n + p 2 p ) + 2A 6 p a Pq . (4) 

p and r are the nuclear density and kinetic energy densityrespectively, and p q and r q are 
the density and the kinetic energy density of species q, respectively. A 1 - A Q are 

Ai = * (1 + ^o), 
A 

A 2 = -to{^ + Xo), 
As = ^[ti(l + -x 1 )+t 2 (l + ^x 2 )}, 
^ = \M^ + x 2 )-t 1 (- + x 1 )], 
A 5 = ^t 3 (l + -x 3 ), 

A 6 = —^t 3 (^ + x 3 ), (5) 
where t , t±, t 2 , t 3 , x , x±, x 2 , x 3 are the parameters of the Skyrme interaction. The effective 



mass m* reads 



2m f 

m* q = m[l + -j^(A 3 p + -A 4 (l ± y)p)]-\ (6) 



Introducing y = Pn , Pp we obtain 

& a pn+pp 



A 3 

Hq(T, P, y) = u q (p, y) + e Coul (p)5 q)P + Tln—p + Tln{\ ± y) 



9s,i 



n+1, , s „, \\ 



6 n (i±y) n (— p)\ (7) 



and 



r 9 = ^ / d 3 kn q (k)k 2 



9. 



s 



I d'kil + expde.-^/T)}- 1 ^ 



= ^ P Y.b n {l±yr\^p)\ (8) 

^ « n=0 

where the symbol " + " stands for neutrons and " — " for protons. Concerning the expansion 

A 3 

in the degree of the degeneracy yj^p) in Pq and r q , in this work, n > 7 term are neglected, 



i.e. one more term than in Refs. [13,14] is included. We find that for temperatures and 
densities relevant to this work(say,T ~ QMeV and p ~ 0.14/m~ 3 ) at least 6 terms are 
needed in order to obtain accurate results. g s is the spin degeneracy (g s = 2) while g s j is 
the spin-isospin degeneracy (g s j = 4). 

Assuming that a hot nucleus is a spherical drop with a sharp edge and uniform distribu- 
tion of nucleons,the Coulomb energy can be expressed as 

QZe 2 3 ,2Zh 2 
£ Coui-^-(^^-. (9) 

Here the Coulomb exchange term is included. Rl is the radius of the nucleus and expressed 
as 

*' = 0- (10) 

The symbols Z and A are the number of protons and nucleons in the liquid drop, respectively. 
The pressure can be obtained through the Gibbs-Duhem relation: 

9P dp n dp p 
Jp- = Pn ^p~ + P *-dp- 



Thus, 



P(T, P ,y) = lj:ji±y)jp^dp : (12) 



q=p,n 

where symbol " + " stands for neutrons and " — " for protons. 

The surface tension term has to be introduced in the pressure for a finite system. For 
the nuclei along the /3-stability line we adopt two different forms of the surface tension. The 
first one is: 

7 (T) = l.UMeV ■ fm'\l + - |)l (13) 

suggested by Ref. [15] and used by Ref. [6,7,9,13,14] in studying the limiting temperature. 
Here we call it Surfl. The other one called Surf2 is: 

7(T) = 7 (0)[(T C 2 -T 2 )/(T C 2 + T 2 )]f (14) 



7 (0) « 18MeV/4irrl r = 1.12/m, (15) 

which was suggested by [16]. This expression of the surface tension is widely used in SMM 
calculations for studying multifragmentation process in heavy ion collisions. It was obtained 
as a parametrization of the calculation of thermodynamic properties of the interface between 
liquid and gas of symmetrical nuclear matter which were performed with Thomas- Fermi and 
H-F method by using Skyrme force. T c is the critical temperature for infinite nuclear matter. 
The pressure caused by the surface tension (Psurf) is given by: 

Psurf = -2 7 (T)/P L . (16) 

Thus, the total pressure of the liquid drop (P) is written as 

P(T,p,y) = P(T,p,y) + P Surf . (17) 

Fig.l shows the isothermal curves of the chemical potential calculated with SLy7 for: 
a)protons with bulk part only, b)neutrons, c)protons with bulk part adding the Coulomb 
term of a hot 208 Pb nucleus. Fig. 2 shows the isothermal curves of the pressure calculated 
with SLy7, in which the Coulomb and the surface tension term are for a hot 208 Pb nu- 
cleus. Fig.3 shows the pressure as a function of density at T=5MeV calculated with SLyl, 
SkM* and SHI in which the Coulomb term and the surface tension term are for a hot 
208 Pb nucleus. There are three sets of curves, namely, the bulk, the bulk+ Coulomb ( 208 Pb) 
and the bulk+Coulomb( 208 P&)+surface tension( 208 P6). In each set, there are three curves 
corresponding to SLy7, SkM* and SHI, crossing at density around p ~ 0.13/m -3 . One 
sees a considerable negative pressure provided by the surface tension through comparing the 
sets with and without surface tension. Furthermore, one can find from the figure that the 
curves for SLyl and SkM* are very close at the densities relevant to this study, while the 
curve for SIII is far away from them except at the crossing point. To show the effect of the 
different forms of the surface tension, in Fig. 4 we plot the pressure at T=5 MeV for 208 Pb 
calculated with the surface tension of Surfl, Surf2, and without a surface tension. One can 
find that the Surfl produces a stronger negative pressure than the Surf2 does, which should 
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have effect on the limiting temperature especially for light nuclear systems. 



III. MASS AND ISOTOPE DEPENDENCE OF LIMITING 

TEMPERATURES 

We adopt the the same model used in refs. [13,14] to study the mass and isotope dependence 
of limiting temperatures for hot nuclei. The model treats the hot nucleus as a uniformly 
charged drop of nuclear liquid with sharp edge at a given temperature, which is in thermal 
equilibrium with the surrounding vapor. The equilibrium between the droplet and the vapor 
surrounding both thermal mechanical and chemical leads to a set of two-phases coexistence 
equation: 

fh(T, p L , y L ) = n P (T, p v , y v ), (18) 
p n (T, p L , y L ) = p n (T, p v , y v ), (19) 

P(T,pL,yL) + Psurf(T,p L ) = P(T, Pv ,yv). (20) 

The subscript L refers to the liquid phase and V to the vapor phase. For simplification, 
the Coulomb interaction in vapor is screened in the calculation of the pressure P{T, p v ,yv) 
and the chemical potential of protons p p (T, pv,yv)- These three coexistence equations with 
three variables can be solved directly to get the coexistence point of the liquid drop and 
the surrounding vapor. By finding the upper boundary of temperature that the coexistence 
equations have solution one obtains the limiting temperature. 

1. Mass Dependence of Limiting Temperatures 

Now let us study the mass dependence of limiting temperatures of nuclei along the /3-stability 
line: 

Z = 0.5A - 0.3 x 10- 2 Al. (21) 
7 



Fig. 5 - Fig. 7 show the mass dependence of limiting temperatures calculated with surface 
tension Surfl and Surf2 and Skyrme interactions SLy7, SkM*, SIII, respectively. The 
data of [10,11] are also plotted in the figures. From these 3 figures one can see that lim- 
iting temperatures decrease with the system size, which is independent of the force used. 
Furthermore, one finds that the limiting temperatures calculated with SLyl and SkM* are 
very close and the largest difference is about or less 0.1 MeV while the limiting temperatures 
calculated with SIII are much higher. This is because the EOS corresponding to SLyl and 
SkM* is of soft one while that for SIII is of stiff one(see Fig. 3). Concerning the surface 
tension, one can see that the limiting temperatures calculated with Surf2 are higher than 
that with Surfl. We find, in general, the calculation results with surface tension Surf2 and 
Skyrme force SLyl agree with experimental data better than other combinations of surface 
tension forms and interaction parameter sets. This agreement may illuminate that the sur- 
face tension may be better described by Surf2 rather than Surfl at T < T c . Therefore, in 
the following calculations we adopt SLyl and Surf2. In Table. 1 we list the values of limiting 
temperatures for the nuclei along the /3-stability line with mass number ranging from 20 to 
250. 

2. Isotope Distribution of Limiting Temperatures 

Fig. 8 shows the isotope distribution of limiting temperatures for nuclei: a)0, b)Ca, and 
c)Zr calculated with Skyrme interaction SLyl and surface tension Surf2, respectively. One 
can see from these figures that the isotope distribution of limiting temperatures appears 
to be a parabolic shape. The centroid of the parabolic curve is not located at the isotope 
on /3-stability line(T=0) but incline to the neutron-rich side. The shift of the centroid of 
the isotope distribution of limiting temperatures to the neutron-rich side seems to be in 
consistent with [17] in which it was mentioned that there was evident that neutron rich 
nuclei had higher limiting temperatures. This effect is due to the fact that the Coulomb 
term does not depend on the temperature explicitly while the isospin dependent part of the 
bulk term does depend on the temperature explicitly, which leads the compensation effect 

8 



of the Coulomb energy and the symmetry energy with respect to the isospin degree to be 
different for hot nuclei and cold nuclei. 

It has been observed that the isospin fractionation distillation effect on multifragmen- 
tation process in neutron-rich heavy ion collisions [18]. We expect that this effect should 
also be shown up in the equilibrium value of the isospin asymmetry of vapor phase at the 
limiting temperature. In Fig. 9 we show the isotope distribution of the equilibrium values 
of the isospin asymmetry of vapor at the limiting temperature, y l ^ m , for O, Ca, and Zr, 
respectively. One can see from the figure that the equilibrium value of y l ^ m increases with 
the increase of the number of neutrons and finally a saturation value might be reached. The 
isospin distillation effect leads vapor to be even neutron-rich for neutron-rich isotopes and 
to be even neutron-lack for neutron-lack isotopes which can be seen in the figure clearly. 
The vanish point of the equilibrium value of y l ^ m , at which the isospin asymmetry of vapor 
changes from negative to positive value, is shifted a little to the neutron-rich side. The 
vanish point should be sensitive to the isospin dependent part of both the interaction and 
the surface tension. 

As is already mentioned in the introduction that the SLy series of Skryme interactions 
is designed to describe the properties of isospin asymmetric nuclei, its isospin dependent 
part is rather different from SkM* [12]. We expect that the isospin distribution of limiting 
temperatures calculated with SLyl should be different from that with SkM* though the 
mass distributions of limiting temperatures calculated with SLyl and SkM* are very close. 
Our calculation results of the isotope distribution of limiting temperature of Zr shown 
in Fig. 10 does demonstrate that the isotope distribution of limiting temperatures depends 
on the isospin dependent part of interaction sensitively. Limiting temperatures calculated 
with SLyl are higher than those with SkM* and the difference between them increases as 
the isospin asymmetry \N — Z\/A increases. It was demonstrated that the SLy series of 
interactions had a better isotope properties than old Skyrme interactions [12], so the isotope 
distribution predicted with SLyl may be more proper. 

As pointed out at Sec. I, the surface tension is a work per surface area needed to create 
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the surface. The symmetry energy term in EOS should also contribute to the surface tension 
and therefore a symmetry-surface tension term should be introduced in the surface tension. 
Here we adopt a surface tension including a symmetry-surface term: 



to study the effect of the isospin dependence of the surface tension. Here a s is taken to 
be 28.5MeV [12,19]. Fig. 11 shows the isotope distribution of limiting temperatures of hot 
nuclei C and O calculated with and without symmetry-surface tension term taken into 
account, respectively. The symmetry-surface tension term raises limiting temperatures of 
neutron-rich isotopes but almost does not affect limiting temperatures at the neutron-lack 
side. Therefore the information of the symmetry-surface tension term can only be obtained 
by the measurement of limiting temperatures of neutron-rich nuclei. 



In summary, in this work we have studied the mass and isotope dependence of limiting tem- 
peratures. The tendency of limiting temperatures decreasing with the increasing of mass 
number for nuclei along the /^-stability line is in agreement with refs. [6,7,9]. The influence 
of different Skyrme forces and different forms of the surface tension on limiting tempera- 
tures of hot nuclei is investigated. We find that the results calculated with newly developed 
Skyrme force SLyl and the surface tension of Surf2 are in good agreement with data. The 
isotope distribution of limiting temperatures appears to be a parabolic shape. The centroid 
of the parabolic curve is not located at the isotope of /3-stability line but inclined to the 
neutron-rich side. The equilibrium value of the isospin asymmetry of vapor phase at the 
limiting temperature y v shows clearly the isospin fractionation distillation effect. An isospin 
dependent surface tension is introduced in studying the isotope distribution of limiting tem- 
peratures. Our study shows that the neutron-rich side of the isotope distribution of limiting 
temperatures depends on the isospin dependence of the the surface tension sensitively. 




(22) 



IV. SUMMARY 
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Figures Caption 



Fig.l The isothermal curves of the chemical potential calculated with SLy7 for: a)protons 
with bulk part only, b)neutrons, c)protons with bulk part adding the Coulomb term of a 
hot 208 Pb nucleus. 

Fig. 2 The isothermal curves of pressure. The bulk part is calculated with SLy7. The 
Coulomb and the surface tension part are of a hot 208 Pb nucleus. 

Fig. 3 The isothermal curve of pressure at T=5MeV calculated with SLy7, SkM*, and 
SIII, respectively. The Coulomb and surface tension term is corresponding to a hot 208 Pb 
nucleus. 

Fig. 4 The comparison of the pressure P at T = 5MeV for three cases, namely, without 
surface tension, with Surfl and Surf2. The surface tension term is corresponding to a hot 
208 Pb nucleus. 

Fig. 5 The mass dependence of the limiting temperatures calculated with Skyrme inter- 
action SLyl and surface tension Surfl and Surf2, respectively. 
Fig.6 The same with Fig.6 but with SKM*. 
Fig. 7 The same with Fig.6 but with SIIL 

Fig. 8 Isotope distributions of the limiting temperatures for O, Ca, and Zr. 

Fig. 9 Isotope distributions of the equilibrium values of the isospin asymmetry of vapor 
at limiting temperature for 0,Ca,Zr, respectively. 

Fig. 10 Isotope distribution of the limiting temperatures for Zr calculated with SLyl and 
SkM*, respectively. 

Fig. 11, the isotope distribution of limiting temperatures of hot nuclei C and O calculated 
with and without symmetry-surface term taken into account, respectively. 

Table. 1 The values of limiting temperatures for nuclei along the /3-stability line. 
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